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Abstract 



m 
o 

Q I I compare the available theoretical models that describe the Cronin effect on 

CN ■ hadron and minijet production in proton-nucleus collisions, pointing out sim- 

ilarities and differences among them. The effect may be summarized by the 
value of two variables. Their values computed in the different models are com- 
pared in the energy range 27.4 GeV - 5500 GeV. Finally, I propose to use the 
pseudorapidity systematics as a further handle to distinguish among the mod- 
els. 



1 INTRODUCTION 
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^ I In this short note I will compare available theoretical models for the description of the so-called Cronin 

^ ' effect [2] in inclusive hadron spectra in proton-nucleus (pA) collisions. The analysis will be limited to 

references containing quantitative predictions in the case of pA collisions [3-7]. The observable which 
I am interested in is the Cronin ratio, R, of the inclusive differential cross-sections for proton scattering 
on two different targets, normalized to the respective atomic numbers A and B: 



o 



A dapB/d?PT 



In absence of nuclear effects one would expect R{pt)=1, but for A>B a suppression is observed experi- 
mentally at small pT, and an enhancement at moderate pT with R{pt) ^ 1 as ^ cxd. This behaviour 
^ \ may be charachterized by the value of three parameters: the transverse momentum px at which R crosses 

' unity and the transverse momentum pM at which R reaches its maximum value Rm=R{pm), see Fig. 1. 

These Cronin parameters will be studied in Sec. 3 

The Cronin effect has received renewed interest after the experimental discovery at RHIC of an 
unexpectedly small < 1 in Au-Au coUisions at y/i=130 GeV compared to pp collisions at the same 
energy [8]. This fact has been proposed as an experimental signature of a large jet quenching suggesting 
that a Quark-Gluon Plasma was created during the collision [9]. However, the extrapolation to RHIC and 
LHC energies of the known Cronin effect at lower energies is haunted by large theoretical uncertainties, 
which may make unreliable any interpretation of signals of this kind. Some light will be shed on this 
problem by the upcoming RHIC data on dA collisions at y^=200 GeV. 

Since in pA collisions no hot and dense medium is expected to be created, a pA run at the same 
nucleon-nucleon energy as in nucleus-nucleus {AA) collisions would be of major importance to test the 
theoretical models and to have reliable baseline spectra for the extraction of novel physical effects. A 
further advantage of pA collisions, is that the multiplicity of particles in the final state is relatively small 
compared to AA collisions. For this reason at ALICE minijets may be observed at transverse momenta 
larger than 5 GeV (see Section on "The experimental parameters for pA at the LHC" in Ref. [1]). As I 
will discuss in Sec. 3, this may allow the use of the Cronin effect on minijet production as a further check 
of the models. 



'Contribution to tlie CERN Yellow report on Hard Probes in Heavy Ion Collisions at the LHC. 



2 THE MODELS 



Soon after the discovery of the Cronin effect [2], it was realized that the observed nuclear enhancement of 
the pt-spectra could be explained in terms of multiple interactions [10, 11]. The models may be classified 
according to the physical object which is undergoing rescatterings 
(the projectile hadron or its partons), and to the "hardness" of the 
rescattering processes taken into account. Note that a parton is com- 
monly said to undergo a hard scattering if the exchanged momentum 
is of the order or greater than approximately 1 GeV. However, phys- 
ically there is no sharp distinction between soft and hard momen- 
tum transfer. Therefore, I prefer to make reference to the so-called 
two component models of hadron transverse spectra, and call hard 
a scattering which is described by a power-law differential cross- 
section at large pt, and soft a scattering whose cross-section is de- 
creasing faster than any inverse power of the transverse momentum 
at large pt. In Tables 2.1 and 2.2 I provide a quick comparison of the hadronic and partonic rescattering 
models, respectively. 




Fig. 1: Definition of Px,Pm,-R A/. 



2.1 Soft hadronic rescattering models [3] [4] 

These models are based on the pQCD coUinearly factorized cross-section for inclusive particle produc- 
tion in pp collisions. In order to describe the large-py tail of transverse momentum spectra, one has to 
include also an intrinsic transverse momentum for the colliding parton. The collision of a proton on a 
nuclear target of atomic number A is then obtained in a Glauber-type model as follows^ : 

^ = kJ2f,/^^ ^.M ® -^(^J - kl) D'l , (1) 

i,j,k,l 

where the proton and nucleus parton distribution functions are, respectively, 

Fi/p = fi/p{xi, Q ) 7T2T and Fj/a = TA{h) fj/p{xj,Q ) . (2) 

In Eq. (1) da/dt{ij kl) is the pQCD parton-parton cross-section, the variables with a hat are the 
Mandelstam variables at parton level, and a sum over incoming and outgoing parton flavours is per- 
formed. The proton is considered point-like compared to the target nucleus, and to scatter on it at impact 
parameter b. The nucleus is described by the Woods-Saxon nuclear thickness function TA{b). In Eq. (2) 
fi/p{A){x, Q"^) are the parton distribution functions of the proton (nucleus); isospin imbalance is taken 
into account and nuclear shadowing is included by the HIJING parametrization [12]. Partons are as- 
sumed to have an intrinsic transverse momentum with average squared value {k'^)pA(Ap) and a Gaussian 
distribution. Due to the ki and kj integrations a regulator mass m=0.8 GeV has been used in the pQCD 
cross-section. Finally, D^{z, Q'^) are the fragmentation functions of a parton k into a hadron h with a 
fraction z of the parton momentum. 

Soft proton-nucleon interactions are assumed to excite the projectile proton's wavefunction, so that 
when the proton interacts with the next target nucleon its partons have a broadened intrinsic momentum. 
Each rescattering of the proton is assumed to contribute to the intrinsic momentum broadening in the 
same way, so that: 

{kT)pA{b, V^) = {k'^)pp + Sx hA{b, V^) , (3) 

'integrations are schematically indicated with crossed circle symbols. For details see the original references. 



Table 2.1 - Parameters of the soft hadronic rescattering models. [4]. 



model 


hard scales GeV 


K 


regul. 


proton intrinsic fer (GeV^) 


n 


average feT-kick (GeV^) 


nPDF 


Soft 


Q = Q' =PT 


2 


0.8 GeV 


(fc2>pp = 1.2 +0.2a4g^)g^ t 


oo 


5f01 0-^^=^ \n^{Q/GeV) 
- U./OO ^^jj^ (Q/GeV) 


HIJING 


Sofl-sal. 




1 


0.8 GcV 




4 


,) = O.i 


HIJING 



^ = 2stu/{s + t + m)^; parametrization chosen to best reproduce pp data. 

^ No explicit parametrization is given. Values of {kT}pp extracted from a "best fit" to pp data, see Fig.15 of Ref. [4]. 
'' These scales used in the computations of Table 3. In Ref. [4] Q =^ and Q' =^ . 



where {k^)pp is the proton intrinsic momentum needed to describe hadron transverse spectra in pp colU- 
sions, 6 is the average momentum squared acquired in each rescattering, and 



[ n 11 VA — i- 



(4) 

> n 



represent the average number of colUsions which are effective in broadening the intrinsic momentum. 
Both models assume Ha to be a function of the number of proton-nucleon collisions Va{^,\fs) = 
app{^/s)TA{b), with app the nucleon-nucleon inelastic cross-section. However, Ref. [3] takes m=(X), 
while Ref. [4] assumes an upper Umit n=4 justified in terms of a proton dissociation mechanism: after 
a given number of interactions the proton is so excited that it can no more interact as a whole with th 
enext nucleon. 1 will call the first model simply soft and the second soft-saturated. In both models target 
nucleon do not have rescatterings, so that 



(4)ap(6,V^) = (4) 



pp 



Further differences between the models are related to the choices of the hard-scales, of the iC-factor 
which simulates NLO contributions to the parton cross-section, and to the parametrizations of {k^)pp 
and 6 (see Table 2.1). 



2.2 Soft partonic rescatterings: the colour dipole model [5] 

In this model the particle production mechanism is controlled by the coherence length lc=y/s/{mNkT), 
where mjv is the nucleon mass and kr the transverse momentum of the parton which fragments in the 
observed hadron. Depending on the value of Ic, three different calculational schemes are considered. 
(a) In fixed target experiments at low energy (e.g., at SPS), where Ic <^ Ra, the projectile's partons 
interact incoherently with target nucleons and high-pt hadrons are assumed to originate mainly from 
projectile's parton which underwent a hard interaction and whose transverse momentum was broadened 
by soft parton rescatterings. The parton is then put on-shell by a single semihard scattering computed 
in factorized pQCD. This scheme I will discuss in detail below, (b) At LHC, where the cm. energy is 
very large and Ic ^ Ra, the partons interact coherently with the target nucleons and high-pr hadrons 
are assumed to originate from radiated gluons; parton scatterings and gluon radiation are computed in 
the light-cone dipole formalism in the target rest frame, (c) At intermediate energies, like at RHIC, an 
interpolation is made between the results of the low- and high- energies regimes discussed above. All the 
phenomenological parameters needed in this model are fixed in reactions different from pA collisions, 
and in this sense the model is said to be parameter-free. 

In the short coherence length scheme, pQCD factorization is assumed to be vahd and formula (1) 
is used with parton masses mg=0.8 GeV and mq=0.2 GeV for, viz., gluons and quarks. Moreover, 

Fi/p = fi/p{xi + ^^,Q^)^ix,b) and F,/^ = r^(6) /../^(x,-, Q^) 6) . 



Table 2.2 - Parmeters of the soft partonic rescattering model (at short-/c) and of the hard partonic 

rescattering models. 



model 


hard scales 


K 


regulators (GeV) 


intr. kr 


dE/dz 


dipole cross-sect. 


nPDF 


Col. dip. 


Q ^ Q' = PT 


® 


mg=0.8, mq=0.2 


as Soft mod. 


-2.5GeV/fm 


(70=23 mb, A=0.288, xo=3-10-'' 


EKS98T 


Hard AT 


Q^Q'^Pt 


2 


/i free param. 


no 


no 


computed from pQCD 


no 


HardGV 


Q^Q' =PT 




^i=0.42^ 


as Soft mod. 


no 




EKS98 



* Q=iji in Ref. [6]. ^ Only at large Xj (EMC effect). ® Factors out in the Cronin ratio. 
/Li determines only the typical momentum transfer in elastic rescatterings. 



Parton rescatterings are computed in terms of the propagation of a qq pair through the target nucleus, and 
the final parton transverse momentum distibution dNi/cPkiT is written as: 



zT 



(dh,d\2e'^^^^^-'^^^ [iM)e-5(-f+-i)(fce>l \e-h<ir,-r^^-)TA(i^)]='^ + 0{agTA) 

J TT L J d^feT 



(5) 



The first bracket in Eq. (5) represents the contribution of the proton intrinsic momentum. The sec- 
ond bracket is the contribution of soft parton rescatterings on target nucleons, expressed through the 
phenomenological dipole cross-secton: for a quark, (r^qirr, x) = ao [l — exp { — Q'l{x))] with 
Qs=i GeV{x/xo)^/'^ (for the value of the parameters see Table 2.2); for a gluon (t|^=9/4 cr|^ is used. 
The expansion of Eq. (5) to zeroth order in gives the intrinsic kx distribution dN^*^^ /d^kx of the 
nucleon partons; the first order term represents the contribution of one-rescattering processes, and so on. 
Partons from the target nucleus are assumed not to undergo rescatterings because of the small size of 
the projectile. Energy loss of the projectile partons is taken into account by a shift of their fractional 
momentum proportional to the energy of the radiated gluons, given by the product of the average mean 
path length AL and the energy loss rate dE/dz [15]. As nuclear shadowing effects are computed theo- 
retically in the dipole formalism, see Eq. (5), parton distribution functions in the target are modified only 
to take into account antishadowing at large x according to the EKS98 parametrization [14]. 

By Fourier transforming the dipole cross-section one sees that the transverse momentum distribu- 
tion of the single parton-nucleon scattering is Gaussian, whence the classification of this model among 
the "soft" ones. However, the single distributions are not just convoluted obtaining a broadening propor- 
tional to the average number of rescatterings. Indeed, in the second bracket the rescattering processes 
have a Poisson probability distribution. As a result, the nuclear broadening of the intrinsic momentum is 
smaller than the product of the average number of rescatterings and the single scattering broadening; this 
might give a dynamical explanation of the assumption used in the Soft-saturated model [4], that n goo 
in Eq. (4). 



2.3 Hard partonic rescattering model [6] [7] 

The model of Ref. [6], hereafter labeled "hard AT", assumes parton rescatterings responsible of the 
Cronin enhancement, and includes up to now in the computations only semihard scatterings, i.e., scatter- 
ings described by the pQCD parton-parton cross-section. It is the generalization to an arbitrary number 
of hard parton rescatterings of the early models of Ref s. [10, 1 1, 16] and of the more recent Refs. [17, 18], 
limited to 1 hard rescattering. As shown in Ref. [11], considering only one rescattering may be a reason- 
able assumption at low energy to describe the gross features of the Cronin effect, but already at RHIC 
energies this might not be enough for the computation of the Cronin peak Rm [6]. The AT model as- 
sumes the S-matrix for a collision of n partons from the on m partons from the target to be factorizable 
in terms of S-matrices for parton-parton elastic-scattering, and assumes generalized pQCD factorization 



[19]. The result is a unitarized cross-section, as discussed in Refs. [6, 20]: 

« i 
The first term accounts for multiple semihard scatterings of proton partons on the nucleus; in the second 
term the nucleus partons are assumed to undergo a single scattering, and = X^i fi\p^'^i\H- Nuclear 
effects are included in dNf^ / cPUt, the average transverse momentum distribution of a proton parton who 
suffered at least one semihard scattering. In impact parameter space it reads 



(7) 



(PkT J 47r 

where unitarity is explicitly implemented at the nuclear level, as discussed in Ref. [20]. In Eq. (7) 
a^^^(r) = KJ2jI d^P [l - e-'P''] f (ij ^ ij) ® Moreover, toi = a^H^Aib) (i.e., the parton- 
nucleon semihard cross-section times the thickness function) is identified with the target opacity to the 
parton propagation. Note that cr^(r) oc as r — >^ and cr^(r) a^jj asr ^ 00. Tliis, together with 
the similarity of Eqs. (7) and (5) suggests the interpretation of cr^(r) as a hard dipole cross-section, 
which accounts for hard parton rescatterings analogously to what tr^ does for soft parton rescatterings. 
Note that no nuclear effects on PDF's are included, but shadowing is partly taken into account by the 
multiple scattering processes. 

To regularize the IR divergences of the pQCD cross-sections a small mass regulator /x is introduced 
in the parton propagators, and is considered a free parameter of the model which signals the scale at 
which pQCD computations break down. As a consequence of the unitarization of the interaction, due 
to the inclusion of rescatterings, both px and pM are almost insensitive on fi [6]. For this reason these 
two quantites are considered a reliable prediction of the model^. Note, however, that they both depend 
on the cm. energy -s/s and on the pseudorapidity 77. On the other hand, Rm is strongly sensitive to 
the IR regulator. This sensivity may be traced back to the inverse-power dependence on ^ of the target 
opacity Ui. oji (xl/ ji"', where the power a > 2 is energy and rapidity dependent. The divergence of Ui as 
ji indicates the need of unitarization of the parton-nucleon cross-section, and deserves further study. 
Therefore, can be here considered only an effective scale which simulates non-perturbative physics not 
considered here [26, 29], the non-linear evolution of parton distribution functions in the target [30] and 
physical effects up to now neglected, e.g., coUisional and radiative energy losses [31]. 

In the model of Ref. [7], hereafter labeled "hard GV", the transverse momentum broadening of a 
parton which undergoes semihard rescatterings is evaluated with the help of Eq. (7) to be 

{kl)H = uJ^?hl{l + c^ , (8) 

where the IR regulator /j, is physically identified with the medium screening mass and fixed to fi = 0.42 
GeV, and represents the typical momentum kick in each elastic rescattering of a hard parton. The factor 
c and the constant term 1 are introduced in order to obtain no broadening for px ^ partons, as required 
by kinematic considerations. The average value in the transverse plane of the opacity uj « (0.4/fm)RA 
(with Ra the nuclear radius) and of the factor c/fi'^ = 0.18 are fixed in order to reproduce the experi- 
mental data at y^=27.4 GeV and y^=38.8 GeV [22]. With these values the logarithmic enhancement 
in Eq. (8) is of order 1 for ~ 3 GeV. Note that uj and c are assumed to be independent of ^/s and r]. 
Finally, the transverse spectrum is computed by using Eq. (1) and adding to the semihard broadening of 
Eq. (8) the intrinsic momentum of the projectile partons, (A;o)=1.8 GeV^: 

with shadowing and antishadowing corrections to target partons as in the EKS98 parametrization [14]. 

^This result is very different from the conclusion of Ref. [18], based on a single-rescattering approximation, that px oc po- 
Hence px cannot be used to "measure" the onset of hard dynamics as proposed in that paper. 



Table 3 - Cronin effect at ri=0: comparison of theoretical models, pm is expressed in GeV. 





model 


charged pions 


partons 


Pm 


Rm 


Ref. 


Pm 


Rm 


Ref. 


21 A GeV 
pW/pBe 
data: Ref. [22] 


Soft 

Soft- saturated 
Color dipole 
Hard AT 
HardGV 


4.0 
4.5*; 4.4© 
4.5 

4 


1.55* 

1.46*; 1.46® 
1.43±0.08® 

1.4 


[3] 
[24] 
[5] 

[7] 


5.1*; 5.1® 
6±0.8t 


1.50*; 1.51® 
1.1^; 1.4'' 


[24] 


200 GeV 
pAa/pp 


Soft 

Soft- saturated 

v^UlUl LllUUlC 

Hard AT 
HardGV 


3.5 
2.9*; 2.7® 
2.7 

3.0 


1.35±0.2t 
1.15*; 1.47® 
1.1 

1.3 


[3] 
[24] 
[5] 

[7] 


4.4*; 4.2® 

7±1^ 


1.29*; 1.70® 

1.25^; 1.2^ 


[24] 


5500 GeV 
pPbjpp 


Soft 

Soft- saturated 
Color dipole 
Hard AT 
Hard GV 


3.5 
2.4*; 2.2® 
2.5 

«40^ 


1.08±0.02f 
0.78*; 1.36® 
1.06 

1.05^ 


[25] 
[24] 
[5] 

[7] 


4.2*; 4.2® 
ll±1.3t 


0.91*; 1.60® 
2.1^; 1.2^ 


[24] 



* With HIJING shadowing [12]. ® Without shadowing. ^ Error estimated by varying dE/dz within error bars [23]. 

* Central value with multiple scattering effects only; error estimated by using different shadowing parametrizations. 
t Numerical errors mainly. ^ Using /j = 1.5 GeV. ^ Using jx = 0.060 -|- 0.283 log(v^), see text. 



^ Completely dominated by EKS98 shadowing and anti-shadowing; result considered not reliable as yet, see text. 
3 PREDICTIONS AND CONCLUSIONS 

In Table 3 I listed the values of the Cronin parameters pm and Rm computed in the various models for 
proton-nucleon center-of-mass energy of 27.4 GeV (representing the low-energy experiments at CERN 
ISR and SPS and at Fermilab), of 200 GeV (RHIC) and of 5500 GeV (LHC). The targets considered in the 
Cronin ratio are listed as well. Since px is of the order of 1 GeV in almost all models at all energies, and 
lies at the border of the validity range of the models, it's value is not shown. Uncertainties of the model 
calculations are included when discussed in the original references (see the notes at the foot of the table). 
In the case of the soft-saturated model, the uncertainty due to the choice of shadowing parametrization is 
illustrated by giving the results obtained with no shadowing, beside the results obtained with the HIJING 
parametrization [12]. Using the "new" HIJING parametrization [13] would change mid-rapidity results 
only at LHC energy, where a 15% smaller Cronin peak would be predicted [24]^. In the case of the "hard 
AT" rescattering model, the major theoretical uncertainty lies in the choice of the effective parameter ^ 
as discussed at the end of Sec. 2.3. In the table, two choices are presented: (a) w energy-independent 
value = 1.5 GeV, which leads to an increasing Cronin effect as energy increases; (b) /lis identified 
with the IR cutoff po discussed in [32], in the context of a leading order pQCD analysis of pp collisions. 
That analysis found po to be an increasing function of ^/s. By performing a simple logarithmic fit to 
the values extracted from data in Ref. [32] we find fi = joo(v^) = 0.060 + 0.283 log (a/s), leading to 
a decreasing Cronin effect. Note that a scale increasing with ^/s appears naturally also in the so-called 
"saturation models" for hadron production in A A colUsions [26-28]. In the "hard GV" model at LHC 
energy, the remnants of Cronin effect at ~ 3 GeV are overwhelmed by shadowing and the calculation 
in this region cannot be considered reliable as yet. The Rm = 1-05 at pT — 40 in Ref. [7] is a result of 
antishadowing in the EK98 parametrization and is not related to multiple initial state scatterings. 

As discussed in the introduction, experimental reconstruction of minijets at ALICE should be 
possible in pA collisions for minijet transverse momenta Pt ^ 5 GeV. The pT-spectrum of the partons 
which will hadronize giving the observed minijet may be obtained by setting D^=5{z — 1) in Eqs. (I) 
and (6). This may be very interesting, because pQCD computations suffer from large uncertainties in 

^Note however that the new parametrization, which predicts a much larger gluon shadowing at a; < 10~^ than the "old" 
one [12], seems ruled out by data on the ratio of Sn and Ca F2 structure functions [34]. 



the determination of FF's at large z, where they are only loosely constrained by existing data [21]. For 
this reason I listed in Table 3 also the Cronin parameters for the case of parton production. However, jet 
reconstruction efficiency should be accurately evaluated to assess the usefulness of this observable. 

Table 3 shows that there are large theoretical uncertainties in the extrapolation of the Cronin effect 
from lower energies to LHC energy. A major source of uncertainty for most of the models is the size of 
nuclear shadowing and anti-shadowing at small x, see Ref. [34] for a detailed discussion and comparison 
of the existing parametrizations. For example the HUING parametrizations [12, 13] predict more gluon 
shadowing than the EKS98 [14] at small x < 10~^. At LHC this is the dominant region at mid-rapidity 
and medium-small transverse momenta. On the other hands the HIJING parametrizations predict less 
antishadowing than EKS98 at x > 10~^, which is the dominant region at large enough pT at all energies. 
At LHC all these effects may lead up to a factor 2 uncertainty in the height of the Cronin peak Rm- 

In conclusion, a pA run at LHC is necessary both to test theoretical models for particle production 
in a cleaner experimental situation than in AA collisions, and to be able to make reliable extrapolations 
to A A colhsions, which is the key to disentangle known effects and new physics. Since, as we have 
seen, the nuclear effects are potentially large, it would be even preferable to have a pA run at the same 
energy as the AA run. In addition, the A systematics, or the study of collision centraUty cuts, would be 
interesting since would allow to change the opacity of the target - then the size of the Cronin effect - in a 
controllable way. Finally, let me remark that the ?7-systematics of the Cronin effect has been considered 
in the literature only in Ref. [6]. However, as discussed also in [33], given the large pseudorapidity 
coverage of CMS this observable might be a very powerful handle for the understanding of the effect. 
It would allow to systematically scan nuclear targets in the low-x region, and would help to test the 
proposed models, in which the rapidity affects the size of the Cronin effect in rather different ways. 
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